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The semiconductor device modeling 
program PC-ID and the programs that 
support its use in high-accuracy model- 
ing of photodiodes, all of which were 
described in Part I of this series of pa- 
pers, are used to simulate the interpo- 
lation of high-accuracy internal 
quantum-efficiency calibrations in the 
spectral region between 450 nm and 
850 nm. Convenient interpolation for- 
mulae that depend only upon wave- 
length arc derived. Uncertainty spectra 
for a number of sources of error are 
also derived. The formulae are normal- 
ized to experimental internal-quantum 
efficiency calibrations in the 440 to 470 



nm spectral region and at 860 nm and 
arc used to interpolate the calibration 
values between these wavelengths. The 
results of the interpolations are com- 
pared with experimental calibration 
data that arc available at a few wave- 
lengths between 440 and 860 nm. The 
disagreement between the interpolated 
and measured internal quantum-effi- 
ciency data is never worse than 0.0003. 
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1. Introduction 

Part II of this series of papers suggested that the 
Hamamatsu^ 1337 type photodiode might be conve- 
niently used for interpolating or extrapolating high- 
accuracy quantum efficiencies in the 400 to 900 nm 
spectral region. The idea of quantum-efficiency ex- 
trapolation was implicit in the first detailed descrip- 
tion of how a silicon photodiode could be used (in 
this case in conjunction with a thermal detector 
having a spectrally flat responsivity) to realize a 
scale of absolute spectral responsivity [1], and inter- 
polation was explicitly considered in reference [2], 
Beyond these two publications, however, little pro- 
gress has been reported in this area until recently. 



' Certain commercial equipment, instruments, or materials are 
identified in this paper to specify adequately the experimental 
procedure. Such identification does not imply recommendation 
or endorsement by the National Institute of Standards and 
Technology, nor does it imply that the materials or equipment 
identified are necessarily the best available for the purpose. 



Hoyt et al. [3] (using early results of the work 
described here) used three Hamamatsu 1337 type 
photodiodes to extrapolate a quantum-efficiency 
calibration obtained at 441.6 nm with a UDT QED 
100 radiometer to 633.0 nm, and reported an aver- 
age difference of 0.04% compared to independent 
measurements based on their electrically calibrated 
cryogenic cavity radiometer. Even though 0.04% is 
quite good by conventional radiometric standards, 
it was a factor of two larger than the combined 
uncertainty estimate for the two measurement 
techniques being compared, and further study was 
proposed by Hoyt et al. More recently, Zalewski 
and Hoyt [4] have reported a more direct test of 
the accuracy of the use of the Hamamatsu 1337 
type photodiode to extrapolate quantum-efficiency 
values. The results of this test agreed to within the 
combined estimated uncertainty, and two errors 
were discovered in the earlier work. When these 
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errors were corrected, agreement within the esti- 
mated uncertainty of the earlier intercomparison 
was obtained. 

The purpose of Part III of this series of papers is 
to develop formulae and uncertainty estimates for 
interpolating and extrapolating internal quantum- 
efficiency calibrations with Hamamatsu 1337 
photodiodes. It will be shown that this type of 
photodiode has characteristics that suit it particu- 
larly well for this task over the 440 to 860 nm spec- 
tral region. The remainder of Part III is organized 
as follows: Section 2 investigates the errors associ- 
ated with the use of relations that are only approxi- 
mate for the 1337 type photodiode, but which 
simplify interpolation and extrapolation of the in- 
ternal quantum deficiency (one minus the internal 
quantum efficiency) for this type of photodiode. 
Section 3 investigates the various sources of error 
associated with the actual calculation of internal 
quantum-deficiency spectra using the semiconduc- 
tor-device modeling program PC-ID and the sup- 
port programs described in Part I. Finally, section 4 
compares the results obtained here with experi- 
mental measurements and demonstrates agree- 
ment within 0.0003. 

2. Errors Associated with the Quantum- 
Deficiency Approximations 

For the purposes of extrapolating and interpolat- 
ing the internal quantum efficiency of Hamamatsu 
1337 type photodiodes, the internal quantum defi- 
ciency can be approximated by a function 5(A,A^ss, 
S,Tt) of wavelength A and of three recombination- 
related variables: Nss the charge number density 
stored in the front-surface oxide, S the surface re- 
combination velocity at the oxide-silicon interface, 
and Ti the minority-carrier lifetime in the rear re- 
gion. 

For silicon photodiodes, S{\,Nss,S,T[) can be 
approximated by a function that has properties that 
make it particularly useful for extrapolating and 
interpolating quantum-deficiency measurements. 
The approximate equation for S(\,Nss,S,Tt) is 



S(\,N,.,S,T,)=K((Nss,S) 8t(\) 



(1) 



+K,(t,) 5,(A), 
where 

5r(A) = 0, (2) 

for A less than some wavelength Aa. 



To see how this is useful for interpolating quan- 
tum-deficiency values, suppose that S(\{,Nss,S ,Tr) 
has been measured at wavelength Af< Ao, and 5(Ar, 
Nss,S,T,) has been measured at wavelength Ar>Ao. 
The values for Ks(Nss,S) and Kc(t,) can then be de- 
termined from 



JiT/xr as. S(\f,Nss,S,Tc) 



and 



^, ^ S(KNss,S,T,)-K!(Nss,S)5fM 



(3a) 



(3b) 



and these values can be used in eq (1) to interpo- 
late the quantum deficiency between Af and A,. In 
the above equations, 5f(A) and 5r(A) are functions 
only of wavelength, Af and A, are the wavelengths at 
which the internal quantum deficiency is to be nor- 
malized, and Ao will generally fall between 600 and 
650 nm. The subscript f stands for the front of the 
photodiode since this is the location of the recom- 
bination mechanisms that cause 3f(A) to exceed 
zero. The subscript r stands for the rear of the 
photodiode since this is the location of the recom- 
bination mechanism that cause 5c(A) to exceed 
zero. Therefore, Sf(\) and 5r(A) will be called the 
front-region and rear-region internal quantum-de- 
ficiency spectra, respectively. The functional de- 
pendence of Ks and Kr will be dropped from the 
notation for the remainder of this paper. 

Exactly how accurate eqs (1-3) are depends 
upon the particular type of photodiode and the 
range covered by the parameters. Figure 1 illus- 
trates the accuracy of eq (1) for the 1337 type 
photodiode (modeled as described in Part II of this 
series of papers) over the 400 to 900 nm spectral 
region for an internal quantum-deficiency spec- 
trum that is a practical upper limit for this type of 
photodiode. Figure 1 compares the internal quan- 
tum-deficiency spectra 5(A,0,71264 cm/s, oo) and 
5 (A ,0,0,1 ms) with their sum, and with the internal 
quantum-deficiency spectrum 5(A,0,71264 cm/s, 1 
ms). The latter spectrum and the sum of the two 
first spectra differ by less than 2x10"^ over the 
entire 400 nm spectral region, and cannot be distin- 
guished at the scale of the figure. (Any errors with 
absolute values less than 2 x 10"^ are negligible for 
the purposes of this paper.) 

Figure 2 illustrates the accuracy of eqs (1-3) over 
the 400 to 900 nm spectral region for the case 
where r, = oo . The spectra in figure 2 were calcu- 
lated for the range of parameters listed in table 1. 
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900 



Figure 1. The quantum-deficiency spectra 5(A , 0,71264 cm/s, 
00 ) (dashed line) and S(A,0,0,1 ms) (dot-dashed line) and their 
sum (solid line). 




- 10 



400 500 600 700 800 
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Figure 2. The quantum-deficiency spectra 5(A,A^ss,5,Tr) for the 
conditions (Case 1 — upper dashed line; Case 2 — lower dashed 
line; Case 3 — upper solid line; Case 4— lower solid line) listed in 
tabic 1. 



Table 1. Recombination related parameters used with the 
doping distributions shown in figure 3 to model the quantum- 
deficiency spectra of Hamamatsu 1337 type photodiodes 



Case 



(cm-') 



S 

(10^ cm/s) 



(s) 



N^ipc) 



1 





0.71264 


00 


dashed curve 


2 


-3X10'2 


1.83500 


00 


dashed curve 


3 


-IX 10" 


100.000 


00 


dashed curve 


4 


-4X10'^ 


1.56000 


00 


soh'd curve 



Cases 1-3 in that table are based on the front- 
region doping distribution Na(x)/Na(0) shown as 
the dashed curve in figure 3, and Case 4 is based 
on that shown as the solid curve in the same figure. 
Of six doping distributions measured for different 
1337 type photodiodes, the two shown in figure 3 
produce the largest change in shape of the internal 
quantum-deficiency spectrum. 




1 2 

DISTANCE (nm) 



Figure 3. The front-region doping distributions Na(x) used to 
calculate the front-region internal quantum-deficiency curves in 
figure 2. 



The shapes of the front-region internal quan- 
tum-deficiency spectra in figure 2 are compared in 
figure 4 by normalizing all of the spectra to 0.01 at 
440 nm and by plotting their differences relative to 
the normalized spectrum for Case 2. Cases 1 and 3 
compared with Case 2 shows how the shape of 
5f(A) changes when its magnitude is increased by 
over a factor of three in association with either a 
substantial decrease or a substantial increase in Ms. 
(Cases 1 and 3 bracket a number of cases that were 
examined in which Nss and 5 were independently 
varied within the range of values shown in table 1.) 

Case 4 compared with Case 2 shows how the 
shape of 5f(A) changes when the dopant distribu- 
tion is changed from the dashed curve to the solid 
curve in figure 4. This effect is much larger than 
that associated with varying Ms and 5. 

The three spectra in figure 4 characterize the un- 
certainties to be associated with the diode-to-diode 
variations in Nss, S, and Na{x) when any one of the 
spectra in that figure is used as Sf(A) with /Tf^O.Ol 
in eqs (1) and (3). These spectra should be multi- 
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plied by Kt/O.Ol if KiP^O.Ol. For the purposes of this 
paper, the 5(A,//ss,5,oo) spectrum corresponding 
to Case 4 in table 1 will be used for St(X) in eqs (1) 
and (3). 




400 500 600 700 800 900 
WAVELENGTH (nm) 

Figure 4. Differences between the internal quantum-deficiency 
spectra of figure 2 for Cases 1 (dashed line), 3 (dotted line), and 
4 (solid line) when normalized to 0.01 at 440 nm, and that for 
Case 2 when normalized to 0.01 at 440 nm. 



The front-region internal quantum deficiency 
also depends upon two front-region recombination 
mechanisms not considered above. These are the 
Auger and band-to-band mechanisms that become 
important in the heavily doped front region of the 
photodiode. Both effects were included in all of the 
spectra described so far by using the default cross 
sections that are built into PC-ID for these recom- 
bination mechanisms. The PC-ID default values 
for the Auger recombination cross sections give an 
internal quantum-deficiency spectrum with a value 
of 9 X 10"* at 400 nm. The shape of this spectrum is 
well modeled by 5f(A) [5]. In fact, an uncertainty of 
±100% can be tolerated in the Auger cross sec- 
tions without causing an error as large as 2 times 
10"^ anywhere within the 400 to 900 nm spectral 
region. 

For the doping profile shown as the solid curve 
in figure 3, the PC-ID default values for the band- 
to-band recombination cross sections give the life- 
times shown in figure 5 as a function of distance 
from the oxide-silicon interface in the photodiode. 
These lifetimes give an internal quantum-defi- 
ciency spectrum with a value of 6 x 10"^ at 400 nm. 
This spectrum is not as well modeled by 5f(A) as is 
the spectrum for Auger recombination, but an un- 



certainty of ±50% in the band-to-band recombina- 
tion cross sections can be tolerated without raising 
the error associated with this uncertainty to 
2x10"^ anywhere in the 400 to 900 nm spectral 
region. 




50 100 150 200 

DISTANCE (nm) 



250 



Figure 5. The default, band-to-band recombination lifetime 
profile calculated by PC-ID for the doping distribution shown 
as the solid curve in figure 3. 

Now consider 5r(A). Figures 6 and 7 illustrate 
the accuracy of eqs (1) and (3) over the 400 to 900 
nm spectral region for the case where 5=0. (Note 
that when 5=0, the internal quantum-deficiency 
spectrum is independent of N^s.) The spectra in fig- 
ure 6 were calculated for t, = 1 ms and 10 ms, with 
Nss=S=0. For the purposes of this paper, S,(X) 
will be set to zero in eqs (2) and (3) for A<610 nm, 
and the spectrum for Tr = l ms will be used for 
5,(A)for A>610nm. 

The shapes of the spectra in figure 6 are com- 
pared in figure 7 by normalizing them to 0.002 at 
860 nm and plotting their differences relative to 
the spectrum for Tr= 1 ms. The normalization value 
of 0.002 was chosen as a practical upper limit. The 
difference spectrum in figure 7 characterizes the 
uncertainties to be associated with the variations in 
rear-region lifetime when the spectrum for Tr = l 
ms is used as 8,(X) with ^r = 0.002 in eqs (2) and 
(3). This spectrum should be multiplied by ii:r/0.002 
ifii:r;^0.002. 

A convenient function of wavelength has been 
fitted to the spectrum corresponding to Case 4 in 
table 1, and another to the spectrum for Tr = 1 ms in 
figure 6. These functions are 



5f(A)=^o [;^(A)-t-^,^^(A)-t-^2y(A)], 



(4) 
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400 500 600 700 800 900 
WAVELENGTH (nm) 

Figure 6. The internal quantum-deficiency spectra 5(A,JVss,'S, Tr) 
for T,= l ms (solid line), and Tr=10 ms (dashed line) with 
JV«=5=0. 




600 650 700 750 800 850 900 
WAVELENGTH (nm) 

Figure 7. Differences between the internal quantum-defi- 
ciency spectra of figure 6 when normalized to 0.002 at 860 nm. 

and 



The appropriate values for Ao, Ai, Ai, Ao, and A2 
for eqs (4) and (5) are listed in table 2. The func- 
tion in eq (5) becomes negative for A < 620 nm, in 
which case eq (5) is replaced with 5r(A) = as men- 
tioned earlier. The values of .^0 in table 2 normal- 
ize Sf(A) to 0.01 at 440 nm, and 5r(A) to 0.002 at 
860 nm, respectively. 

With Af= 0.5516405, eq (4) fits the spectrum for 
Ms= - 3 X 10^^ cm/s in figure 4 with a residual stan- 
dard deviation of 2.9 x 10"* over the 420 to 900 nm 
spectral region and never differs from that spec- 
trum by more than 7 x 10~* over the same spectral 
region. This translates to less than 10 ppm when 
the internal quantum deficiency is normalized to 
0.01 at 440 nm; thus the error associated with the 
use of eq (4) is negligible for the purposes of this 
paper. The larger variations below 420 nm are neg- 
ligible with respect to uncertainties already identi- 
fied and others discussed below. 

With ^r= 1.685081, eq (5) fits the spectrum 
shown as a solid line in figure 6 with a residual 
standard deviation of 1.4 x 10"^ over the 610 to 900 
nm spectral region, and never differs from that 
spectrum by more than 2.5 x 10"^ over that same 
spectral range. This translates to less than 
1.7 X 10~^ when the quantum deficiency is normal- 
ized to 0.002 at 440 nm; thus the error associated 
with the use of eq (5) is also negligible for the pur- 
poses of this paper. 



Table 2. Values of the parameters to be used in eq (4) for 
extrapolating a short wavelength internal quantum-deficiency 
value to longer wavelengths and in eq (5) for extrapolating a 
long wavelength internal quantum-deficiency value to shorter 
wavelengths 



Parameter 


eq(4) 


eq(5) 


Aa 


1.38002 xlO^ 


2.02157 xl0-« 


Ai 


1.47530x10^ 


3.51764x10' 


At 


8.07476x10-5 


0.0 


As 


4.39312x10' nm 


7.53991 X 10' nm 


A2 


4.27998x10^ nm 


1.66889 X 10" nm 



8r{\)=Aa[llX{\)~A^{k)], 
respectively, where 

;^(A) = exp(-A/Ao), 
and 

y(A) = exp(-(A/A2)^). 



(5) 



(6) 



(7) 



3. Errors Associated with 
Deficiency Calculations 



Quantum- 



Two types of error associated with the quantum- 
deficiency values calculated by PC-ID can be 
distinguished: 1) numerical errors associated with 
the algorithms used by PC-ID and the limited 
number of finite elements available to PC-ID, and 
2) physical errors associated with approximations 
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and simplifications in the physical models used 
with, or built into, PC-ID. Both types are consid- 
ered in this section. 

It was shown in reference [5] that 5f(A) is well 
approximated by 



.(.._ r{I{0)-I[a{X)]} 
^(^)- i + ri[a(X)] ' 

/[a;(A)] k(x) exp(-a;(A>c)dc, 

Jo 

k(x)=k(x)/k(0), 
k(x) = l/[mo(x)D„,(x)l 
r=S/D„(0), 



(8) 

(9) 

(10) 
(11) 
(12) 



where a (A) is the absorption-coefficient spectrum 
of silicon, x is the distance from the oxide-silicon 
interface toward the rear of the photodiode, mo(x) 
is the equilibrium minority-carrier concentration at 
the point x in the front region of the photodiode, 
and D„,(x) is the equilibrium diffusion constant for 
the minority carriers at the point x. Note that the 
sign in the numerator in eq (8) is correct and that 
the sign is incorrect in reference [5]. Also, note 
that /(A) becomes very insensitive to xi as xi in- 
creases beyond some critical value. The point 
where the electron and hole concentrations cross is 
a convenient choice for^ri, but smaller values would 
work as well. 

From the point of view of PC-ID, the kernel 
k{x) in eq (11) is given by 



k(x) =qMo(x)/lnii(x)kTiM!a(x)], 



(13) 



where Mo(x) is the equilibrium majority carrier 
concentration in the front region, rtieix) is the ef- 
fective intrinsic carrier concentration [6], which in 
PC-ID includes the effect of Fermi-Dirac statistics, 
iim(x) is the minority-carrier mobility as a function 
of position in the front region, k is the Boltzmann 
constant, and T is the temperature of the photodi- 
ode. 

Figure 8 compares exp(-a;(400 nm)r) with the 
equilibrium functions Mo(x), nic(x), and (JLm(:«:) as 
calculated by PC-ID over the range 0< a: <200 nm 
for the same photodiode model used to calculate 
the quantum-deficiency spectra in figure 3 for Case 

2 of table 1. Equations (8-12), figure 8, and figure 

3 show why the shape of 5f(A) is very insensitive to 
Nss while being quite sensitive to the shape of 
Na(x)/Na(0). Since the width of the accumulation 
layer created by A''ss is small compared to the dis- 



tance over which exp(-a;(A)x) changes signifi- 
cantly for A >400 nm, changes in A^ss have the effect 
of multiplying the integral in eq (9) by a constant 
factor independent of A>400 nm. On the other 
hand, the changes in Na(x)/Na(0) shown in figure 3 
extend well into the region where exp(-a(\)x) 
varies significantly from unity for A = 400 nm. As A 
is increased above 400 nm, the shape of 5f(A) will 
become less dependent on the shape of Na(x)/ 
Na(0) because a(\) decreases with increasing A. 
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Figure 8. Comparison of itifoC*) (dotted line), h\e{x) (dashed 
line), )x (dot-dashed line), and exp(-a(400 nm) j:) (solid line) 
for a 1337 type photodiode with N^, S, and r, set as defined in 
Case 2 of table 1. The tildes in the preceding expressions indi- 
cate that these expressions have been normalized to unity at 
j: =0, as shown for k(x) in eq (10). 

The numerical accuracy with which PC-ID cal- 
culates Mo(x) within the framework of the drift- 
diffusion approximations for the conditions of 
charge accumulation at an oxide-silicon interface 
as shown in figure 8 has been studied previously 
and found to he quite good [7]. With a uniform 
doping density of 10" cm~^ and an oxide charge 
density of 10" cm~^ PC-ID overestimates Mq(x) by 
about 9% atA: = 0, decreasing to within ±1% for 
x>l nm. The uncertainty ink(x) due to this source 
of error will be modeled by replacing k(x) by 



ki(x)=k(x)/(l-0.1x/xo), 

for X <xo, and 
ki(x)^k(x)/0.9, 



(14) 



(15) 



486 



Volume 96, Number 4, July-August 1991 

Journal of Research of the National Institute of Standards and Technology 



for X >xo in the kernel of the integral in eq (9), 
where xo = 1 nm. 

Another source of error that must be considered 
is the numerical accuracy with which PC-ID would 
calculate the steady-state photocurrent if the equi- 
librium carrier concentration were exact. The un- 
certainty associated with this source of error was 
estimated as follows: 1) The internal quantum- 
deficiency spectrum was calculated directly from 
the PC-ID solution of the steady-state drift- 
diffusion equations for the total current flowing in 
the photodiode for the conditions described by 
Case 2 of table 1. 2) Equations (8-13) were used 
to calculate the internal quantum-deficiency 
spectrum from the equilibrium values of Mo(x), 
ttieix), and yLm(^) calculated by PC-ID for the same 
conditions as in 1) above. The integration in eq (9) 
was carried out using a generalization of Simpson's 
rule for nonequidistant points. 3) Both spectra 
were normalized to 0.01 at 440 nm, and the differ- 
ences calculated. There is no reason to believe that 
these differences, which exceed 2 x 10"^ only at 400 
and at 900 nm, are caused by errors associated with 
PC-ID. It is just as likely that they are associated 
with eqs (8-12) since these equations were derived 
from a number of assumptions that are not rigor- 
ously satisfied. The important point is that it is un- 
likely that either calculation is in error by much 
more than the differences between them. This veri- 
fies the numerical accuracy of the photocurrents 
calculated by PC-ID to the level of accuracy re- 
quired in this paper. 

There are a number of other errors, which are 
associated with incomplete or approximate physical 
models built into PC-ID, that also affect the accu- 
racy with which PC-ID calculates the functions 
defining k{x) in eq (13). These include quantum- 
mechanical effects (tunneling) [8] that change the 
shape of Mo(jc) from that calculated from the drift- 
diffusion approximation for x<\ nm, band-gap 
narrowing effects [9] associated with heavy doping 
and large surface-fields [10] that cause n-Kix) [6] to 
depend upon position in the photodiode, and car- 
rier-carrier and carrier-ion scattering mechanisms 
that cause p^ix) to depend upon position [11]. It is 
beyond the scope of this paper to investigate each 
of these sources of error in detail, but the uncer- 
tainty in k{x) due to this source of error will be 
modeled by replacing k{^) by 



A:2(jc)=A:(jc){l-0.1 log[^(jc)]} 



(16) 



in the kernel of the integral in eq (9). This allows 
an error of 10% per decade change in k{x). 



The differences between the internal quantum- 
efficiency spectrum calculated from eqs (8-13) by 
replacing A:(x) with either ki{^) or fo(x) in eq (9) 
and the spectrum calculated using k{x) were calcu- 
lated after all three spectra were normalized to 
0.01 at 440 nm. Neither difference exceeded 
2 X 10"^ over the entire 400 to 900 nm spectral re- 
gion. Therefore, all of these sources of error are 
considered negligible for the purposes of this pa- 
per. 

Another source of error that is obvious in eqs 
(8-12) is the uncertainty associated with the ab- 
sorption-coefficient data used in the simulations. 
As described in Part I, the absorption-coefficient 
data were calculated from an equation [12] fitted 
to the data of Weakliem and Redfield [13]. An al- 
ternate set of absorption-coefficient data was de- 
scribed by Philipp [14]. Figure 9 plots the 
differences between the internal quantum-defi- 
ciency spectra normalized to 0.01 at 440 nm for 
Case 2 of table 1 when based on the absorption-co- 
efficient data in reference [14] and when based on 
the equation in reference [12]. These differences 
will be used as the estimated uncertainties arising 
from the uncertainties in the silicon absorption-co- 
efficient spectrum when 5f(A) is used to extrapo- 
late an internal quantum deficiency of 0.01 from 
440 nm to longer wavelengths with a 1337 type 
photodiode. A similar calculation was carried out 
for the case where 5=0 and r, = 1 ms. Both spectra 
were normalized to 0.002 at 860 nm. The differ- 
ences, which are shown in figure 10, are taken as 
the estimated uncertainties arising from the uncer- 
tainties in the silicon absorption-coefficient spec- 
trum where 5r(A) is used to extrapolate an internal 
quantum deficiency of 0.002 from 860 nm to 
shorter wavelengths with a 1337 type photodiode. 

Figures 11 and 12 plot the quadrature sum of the 
difference spectra plotted in figures 4 and 9 and in 
figures 7 and 10, respectively. The spectrum in fig- 
ure 11 plots the estimated (one standard deviation) 
uncertainty as a function of wavelength that is asso- 
ciated with the use of 5f(A) in eq (4) to extrapolate 
a measured internal quantum deficiency of 0.01 
from 440 nm to longer wavelengths. The major 
sources of error contributing to the uncertainty 
spectrum in figure 11 are the diode-to-diode varia- 
tions in Na{x)/Na(0) and the uncertainty in the ab- 
sorption-coefficient spectrum for silicon, but a 
number of other sources of error were identified 
with the help of eqs (8-12), and any that produced 
differences greater than 2x10"' for 400 nm 
< A <900 nm are included in figure 11. 
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500 600 700 800 
WAVELENGTH (nm) 



900 



Figure 9. Differences between the internal quantum-deficiency 
spectrum for a 1337 type photodiode satisfying the conditions de- 
fined in table 1 for the absorption-coefficient data of reference 
[14] and for the absorption-coefficient data of reference [12]. 
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Figure 11. Sum in quadrature of the curves in figures 4 and 9. 
This curve is the nominal uncertainty associated with the use of 
eq (4) in eq (2) to extrapolate a quantum deficiency of 0.01 at 440 
nm to any other wavelength between 400 and 900 nm. Notice 
that the sum is plotted on different scales below and above 440 
nm. 
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Figure 10. Differences between the internal quantum-deficiency 
spectra for a 1337 type photodiode with 5=0 and Tr=l ms for 
the absorption-coefficient data of reference [14] and for the ab- 
sorption-coefficient data of reference [12]. 
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Figure 12. Sum in quadrature of the curves in figures 7 and 10. 
This curve is the nominal uncertainty associated with the use of 
eq (5) in eq (2) to extrapolate a quantum deficiency of 0.02 at 860 
nm to any other wavelength between 400 and 900 nm. Notice 
that the sum is plotted on different scales below and above 440 
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The spectrum in figure 12 plots the estimated 
(one standard deviation) uncertainty as a function 
of wavelength that is associated with the use of 
5r(A) in eq (5) to extrapolate a measured internal 
quantum deficiency of 0.002 from 860 nm to 
shorter wavelengths. The major sources of error 
contributing to the uncertainty spectrum in figure 
12 are the change in shape of Sr(A) with xr and the 
uncertainty in the absorption-coefficient spectrum 
for silicon. All other sources of error considered 
were negligible with respect to 2x 10"' for 400 nm 
<X.<900nm. 

No equivalents to eqs (8-12) were used to guide 
the error analysis for 5r(A). The results of the un- 
certainty analysis of 5f(A) show that numerical ac- 
curacy of PC-ID will be a negligible source of 
error. This leaves the errors associated with the ap- 
proximations in the physical models that are built 
into PC-ID. Some of these were tested, but no 
tests were possible for the use of a single SRH trap 
level to model the recombination in the rear region 
of the photodiode. This approximation could be a 
non-negligible source of error; it was shown in Part 
I that this approximation introduced errors of the 
order of the effect being modeled in nonlinearity 
simulations on a UV444B type photodiode. A more 
complete model might make t, depend upon posi- 
tion in the photodiode, which might modify the 
shape of 5r(A) and add an uncertainty that is not 
included in figure 12. 

4. Comparison with Experiment 

To use the results developed thus far, it is neces- 
sary to define the internal quantum-deficiency in- 
terpolating function 5(A), and to assume that it 
exactly satisfies eqs (1-3), which can be rewritten 
in more compact form as 



5(A) = 



&(Af) 



¥A) 
5f(Af) 



+ 



,(4_|ga]4(Ai 



5r(A0 J d,i\,y 



(17) 



where 6i(Af) and 5x(Ar) are the measured values of 
the internal quantum deficiency at the wavelengths 
Af and Ar, respectively; 5f(A) and 5r(A) are defined 
in eqs (4) and (5), respectively, and the constants 
used in these equations are defined in table 2. The 
uncertainty associated with 5(A) due to the sources 
of error considered in the last two sections of this 
paper can be obtained by summing in quadrature 
each of the terms on the right-hand side of the dif- 
ferential of eq (17). That differential is given by 






+ 



"HWI] 



^"Hmib-m 



(18) 



The differentials dSt(\) and ddc(\) are plotted in 
figures 11 and 12 for Af=440 nm, and for Ar = 860 
nm, respectively. 

There are not many data available against which 
to test eqs (17) and (18), but there are some. For 
instance, Zalewski and Hoyt [4] have reported the 
internal quantum efficiencies at 441.6 nm and the 
spectral responsivities and absorptances at 633.0 
nm for five multiple reflection (trap) [4,15-17] ra- 
diometers based on Hamamatsu 1337 photodiodes 
at 441.6 and 633.0 nm. Similarly, Fox [17] has re- 
ported the average internal quantum efficiency of 
10 trap radiometers also based on Hamamatsu 
1337 photodiodes at six Ar"^ ion laser lines, and he 
has also measured the oxide-bias correction for a 
single Hamamatsu 1337 photodiode at four Ar"^ 
ion laser lines. Zalewski and Hoyt [4] presented 
their data, which are summarized in bold face type 
in table 3, as an intercomparison of spectral re- 
sponsivity. Since the internal quantum deficiency is 
the quantity of interest in this paper, it is necessary 
to compute this quantity from the data reported by 
Zalewski and Hoyt as shown in table 3. The correc- 
tions for nonlinearity that Zalewski and Hoyt 
derived from reverse-bias measurements also elimi- 
nate any effects of recombination in the rear region 
of the photodiodes. 

Table 3. Average values of, and uncertainty estimates for, the 
internal quantum efficiency Cf, spectral responsivity R, one mi- 
nus the reflectance (1 — p) and nonlinearity correction (1 - NL) 
reported in reference [4] (bold face type) at wavelength A for five 
multiple-reflection radiometers based on 1337 type photodiodes, 
and internal quantum-deficiency values d,(A) calculated from 
the reported data 



A 
(nm) 



R 
(A/W) 



(1-p) (1-NL) S,(A) 



441.6 +0.99347 
±0.00040 



+ 0.00653 
±0.00040 



633.0 +0.50840 +0.99704 +0.99972 +0.00097 

±0.00009 ±0.00014 ±0.00004 ±0.00023 
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Fox's [17] oxide-bias data, which are summarized 
in table 4, were measured at 10 V, and were re- 
ported as a fractional increase in photocurrent / 
given by •yo=/(10V)//(0) — 1. The oxide-bias exper- 
iment reported in Part II shows that the internal 
quantum efficiency at zero bias Sx is given by 



& = 1-1/(1 + 1.13 To), 



(19) 



where the maximum oxide bias-voltage is 10 V. The 
average internal quantum efficiencies reported by 
Fox are also shown in table 4, and the internal 
quantum deficiencies calculated from them are 
also shown there. 



mated uncertainty significantly because the terms 
are added in quadrature, and the first term on the 
right-hand side of eq (18) is larger than the second 
term for both sets of data. Fox's data point at 
406.85 nm was included in table 6 even though the 
uncertainties associated with extrapolating Si(\) to 
wavelengths shorter than Af are very unfavorable 
for high-accuracy applications. The good agree- 
ment at this wavelength must be considered fortu- 
itous considering the uncertainty associated with 
the predicted value. The differences at the other 
wavelengths never exceed 0.00015, and they fall 
within the estimated uncertainty for the differ- 
ences. 



Table 4. Oxide-bias data ti, reported in reference [17] for a 
single 1337 type photodiode, the average internal quantum-effi- 
eiency data Ct reported in reference [17] for ten multiple-reflec- 
tion radiometers based on 1337 type photodiodes, and the 
quantum deficiencies 5, (A) calculated from these sets of data 



A 
(nm) 



oxide bias 
ro S,(A) 



iqe 



S.(A) 



406.85 


+ 0.0126 


+ 0.01404 


+ 0.9867 


+ 0.0133 




±0.00015 


±0.00017 


±0.0002 


±0.0002 


468.18 


+ 0.0044 


+ 0.00495 


+ 0.9934 


+ 0.0066 




±0.00015 


±0.00017 


±0.0003 


±0.0003 


568.35 






+ 0.9977 
±0.00020 


+ 0.0023 
±0.0002 


647.30 


+ 0.0007 


+ 0.00079 


+ 0.9986 


+ 0.0014 




±0.00015 


±0.00017 


±0.0002 


±0.0002 


799.54 






+ 0.9989 
±0.0002 


+ 0.0011 
±0.0002 


859.07 


+ 0.0000 


+ 0.00000 


+ 0.9987 


+ 0.0013 




±0.00015 


±0.00017 


±0.0002 


±0.0002 



Equation (17) was normalized to the results of 
Zalewski and Hoyt at Af= 441.6 nm, and the predic- 
tions of that equation are compared with the re- 
sults of Zalewski and Hoyt at 633 nm in table 5. 
Note that S^(K) = for this data set due to the ap- 
plication of reverse bias. Equation (17) was also 
normalized to the oxide-bias results of Fox at 
Af =468.18 nm, and the predictions of that equation 
are compared with the oxide-bias results in table 6. 
Again note that Sx(Ar) = because oxide-bias mea- 
surements are not sensitive to recombination in the 
rear of the photodiode. 

The uncertainty estimates for eq (17) in tables 5 
and 6 made use of the data in figure 11 even 
though those data apply to Af=440 nm rather than 
441.6 or 468.2 nm. This does not distort the esti- 



Table 5. Comparison of predictions of eq (17) with measured 
internal quantum-deficiency data S,(A) in table 3 for Af=441.6 
nm and 5r(Ar) = 



Prediction 

wavelength 

(nm) 



Measured 

value 
reference [4] 



Predicted 
value 
eq (17) 



Difference 



633.0 



+ 0.00097 
±0.00023 



+ 0.00089 
±0.00007 



-0.00008 
: 0.00024 



Table 6. Comparison of predictions of eq (17) with internal 
quantum-deficiency data 5t(A) derived from oxide-bias data in 
table 4 for Af=468.18 nm and 5,(Ar) = 



Prediction 

wavelength 

(nm) 



Measured 

value 
reference [15] 



Predicted 
value 
eq (17) 



Difference 



406.85 



647.30 



859.07 



+ 0.01404 


+ 0.01321 


-0.00083 


±0.00017 


±0.01144 


±0.01144 


+ 0.00079 


+ 0.00090 


+ 0.00010 


±0.00017 


±0.00030 


±0.00034 


+ 0.00000 


+ 0.00015 


+0.00015 


±0.00017 


±0.00005 


±0.00018 



Table 7 and figure 13 compare the internal 
quantum-deficiency results of Fox [17] with the 
predictions of eq (17) when normalized at 
Af= 468.18 nm and at A, = 859.07 nm. The data in 
figures 11 and 12 that apply for Af =440 nm and for 
Ar = 860 nm, respectively, are used with eq (18) to 
estimate the uncertainties assigned to the predic- 
tions of eq (17). The difference at 406.85 nm is 
much larger than that obtained with the oxide-bias 
data but still falls within the estimated uncertainty. 
Because the oxide-bias experiment is not sensitive 
to the quantum yield for electron-hole pair produc- 
tion, it is tempting to imagine that the discrepancy 
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Table 7. Comparison of predictions of eq (17) 
internal quantum-deficiency data in table 4 for 
and Ar = 859.07 nm 



Prediction 

wavelength 

(nm) 



Measured 

value 
reference [15] 



Predicted 
value 
eq (17) 



with measured 
Af=468.18 nm 



Difference 



406.85 nm 


+ 0.0133 


+ 0.01744 


+ 0.00414 




±0.00025 


±0.01000 


±0.01000 


568.35 nm 


+ 0.0023 


+ 0.00224 


-0.00006 




±0.00030 


±0.00007 


±0.00031 


647.30 nm 


+ 0.0014 


+ 0.00121 


-0.00019 




±0.00020 


±0.00004 


±0.00020 


799.54 nm 


+ 0.0011 


+ 0.00082 


-0.00028 




±0.00021 


±0.00009 


±0.00023 
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Figure 13. Comparison of measured data of reference [17] 
(open circles) with eq (17) (solid line) normalized to pass 
through the measured point at 468.18 and 859.1 nm. At these 
wavelengths the error bars are the uncertainties in the mea- 
sured values; at the other wavelengths the error bars are the 
quadrature sum of the uncertainties in the measured values and 
those associated with the use of eq (17). 

between the oxide bias and internal quantum effi- 
ciency reflects the fact that the quantum yield for 
the 1337 type photodiode is greater than unity at 
406 nm. However, the fact that this discrepancy is 
not statistically significant at the one-sigma level 
shows that this conclusion cannot be drawn from 
the data presented here. The differences at the 
other wavelengths never exceed 0.0003, and only 
fall outside the estimated uncertainties for the dif- 
ferences at 799.54 nm. 



The results shown in tables 5-7 and figure 13 
verify that the internal quantum-deficiency spec- 
trum of multiple-reflection radiometers based on 
Hamamatsu 1337 photodiodes can be interpolated 
over the 440 to 860 nm spectral region from only 
two measured quantum deficiencies, one at each 
end of the region. This greatly reduces the number 
of measurements needed for high-accuracy spectral 
calibrations of these photodiodes and radiometers. 
Furthermore, the theoretical uncertainty analysis 
suggests that the uncertainties in the interpolated 
values will be less than or equal to the uncertain- 
ties in the measured values, at least as long as the 
latter are greater than ±0.0001. The results in ta- 
bles 5 to 7 and figure 13 do not contradict this idea, 
but the experimental data are not accurate enough 
to confirm it. On the other hand, the experimental 
data do confirm that an uncertainty of ±0.0003 is 
obtainable. 



4. Conclusion 

Part III of this series of papers has shown that 
Version 2 of the semiconductor-device modeling 
program PC-ID can be used to model the spectral 
shape of the internal quantum deficiency of 1337 
type photodiodes with very small uncertainty over 
the 440 nm spectral region. The largest uncertain- 
ties are caused by uncertainties in the absorption- 
coefficient spectrum of silicon, diode-to-diode 
variations in the shape of the front-region dopant 
distribution, and diode-to-diode variations in the 
position-independent rear-region minority-carrier 
lifetime. It was not possible to estimate the uncer- 
tainty caused by approximating the rear-region 
minority-carrier lifetime as being independent of 
position. 

Simple formulae were derived that allow high- 
accuracy internal quantum-efficiency calibrations 
on 1337 type photodiodes to be interpolated from 
measured values at the ends of the 440 to 860 nm 
spectral region over the interior of that region. The 
uncertainties that can be obtained with these for- 
mulae are comparable to those that can be 
obtained with the highest accuracy measurements 
currently available. 
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